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We analyse the production of multileptons in the simplest supergravity model with 
bilinear violation of R parity at the Fermilab Tevatron. Despite the small i?-parity 
violating couplings needed to generate the neutrino masses indicated by current 
atmospheric neutrino data, the lightest supersymmetric particle is unstable and can 
decay inside the detector. This leads to a phenomenology quite distinct from that 
of the i?-parity conserving scenario. We quantify by how much the supersymmetric 
multilepton signals differ from the i?-parity conserving expectations, displaying our 
results in the tuq (8) ?t^i/2 plane. We show that the presence of bilinear i?-parity 
violating interactions enhances the supersymmetric multilepton signals over most of 
the parameter space, specially at moderate and large mo. 
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I. INTRODUCTION 

Supersymmetry (SUSY) provides a promising candidate for physics beyond the Standard 
Model (SM). The search for supersymmetric partners of the SM particles constitutes an 
important item in the agenda of current high energy colliders like the Tevatron, and future 
colliders like the CERN Large Hadron Collider or a linear e~^e~ collider. So far most of 
the effort in searching for supersymmetric signatures has been confined to the framework of 
i?-parity conserving realizations Q]; see, however, Ref. |2] and references therein. 

Recent data on solar and atmospheric neutrinos give a robust evidence for neutrino 

n 

conversions 3], probably the most profound discovery in particle physics in the recent years. 

It has been suggested long ago that neutrino masses and supersymmetry may be deeply 

tied together Ij. Indeed, SUSY models exhibiting /2-parity violation can lead to neutrino 
n 

masses and mixings IS] in agreement with the current solar and atmospheric neutrino data. 

riQ 

Furthermore, the simplest possibility is bilinear i?-parity violation [y, [TJ which may arise 
as an effective description of a spontaneous i?-parity violation from some 

suitable ab initio symmetries 

It is interesting to notice that neutrino mass models based on i?-parity violation can be 



tested at colliders 
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|. In this work, we study the production of multileptons (> 3i 
with £ = e or /i) at the Fermilab Tevatron within the framework of the simplest supergravity 
(SUGRA) model without -R-parity conservation [tI]. We show that the presence of bilinear 
i?-parity violating (BRpV) interactions enhances the supersymmetric multilepton signals 
over most of the parameter space, specially at moderate and large mo. In order to make 
the comparison with the i?-parity conserving case easier, we performed our detailed study 
of the signal and SM backgrounds, adopting the same cuts (soft cuts SC2) proposed in Ref. 

Q 

In SUGRA with universal soft breaking terms at unification, the masses of the sleptons, 
the lighter chargino (xf), and the lighter neutralinos {xi ^-nd X2) considerably smaller 
than the gluino and squark ones over a large range of the parameter space [l3^. Therefore, 
the production of charginos and neutralinos provides the largest reach at the Tevatron. 
In i?-parity conserving scenarios, a promising signal for SUGRA at the Tevatron is the 
production of X2xf pahs and their subsequent decay into three charged leptons in association 
with missing energy due to the undetected lightest supersymmetric particle (LSP), which 
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turns out to be Xv the presence of i?-parity violation, the LSP is no longer stable, 
giving rise to events containing more particles, and it can be electrically charged since it 

n n 

decays. Interesting signals for such scenarios have been worked out for staus 12], stops jl4| . 
gluinos and also for supersymmtric Higgs bosons . The larger activity in the detector 
can either enhance the signal, via the production of additional leptons or suppress it due to 
the hadronic activity that spoils the isolation of the leptons. Therefore, a careful reanalysis 
of the trilepton signal is necessary. 

In this work, we consider a SUGRA model that includes the following bilinear terms in 

W^BRpV = W^MSSM — ^ab^iLIHt ) (l) 

where the last term violates R parity. In order to reproduce the values of neutrino masses 
indicated by current ncutnno data 1181 we ,„„st have 1.1 « 1,|. where , denotes the SUSY 

bilinear mass parameter p| . The relevant bilinear terms in the soft supersymmetry breaking 
sector are 

where the terms proportional to Bi are the ones that violate R parity. The explicit -R-parity 
violating terms induce vacuum expectation values (vev) t>j, i = 1,2,3 for the sneutrinos, in 
addition to the two Higgs vev's f„ and Vd- 

Our goal is to determine the impact of -R-parity violation on SUSY multilepton signals, 
for example in the production of three or more electrons or muons, taking into account the 
magnitude of the mass scale indicated by the current atmospheric neutrino experiments. In 
phenomenological studies where the details of the neutrino sector are not relevant, it has 



been proven very useful to work in t 



le approximation where R parity and lepton number are 



violated in only one generation jl6|,|l9|. Thus, the first step to achieve our goal is to assume 
the approximation where R parity is violated only in the third generation (of course all gauge 
interactions are treated in the full three-generation scheme) . This is the theoretically natural 
choice to make, since the third generation forms the basis for the radiative breaking of the 
electroweak symmetry, driven by the top quark Yukawa coupling. From the point of view 
of the analysis presented below, the breaking of R parity only in the third generation also 
corresponds to the worst-case-scenario: the breaking of R parity in the muon channel would 
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produce muons directly, not just as tau decay products, leading to an enhanced multilepton 
(multi-muon) signal. 

II. MAIN FEATURES OF OUR BRPV MODEL 

The parameter space of our SUGRA model, which exhibits -R-parity violation only in the 
third generation, is 

mo , mi/2 , tan/5 , sign(/i) , , , and m,,3 , (3) 

where mi/2 and mo characterise the common gaugino mass and scalar soft SUSY breaking 
masses at the unification scale, Aq is the common trilinear term, and tan/5 is the ratio 
between the Higgs field vev's. Although many parameterisations are possible, it is convenient 
to characterise the BRpV sector by the bilinear superpotential term and the neutrino mass 

We considered the running of the masses and couplings to the electroweak scale, assumed 
to be the top mass, using the one-loop renormalization group equations. In the evaluation of 
the gaugino masses, we included the next-to-leading order (NLO) corrections coming from 
as, the two-loop top Yukawa contributions to the beta-functions, and threshold corrections 
enhanced by large logarithms; for details see j^]. The NLO corrections are especially 
important for M2, leading to a change in the wino mass up to 20%. We then input the soft 
terms into PYTHIA that was used to evaluated the masses and decay rates of all particles, 
except the first and second neutralinos. 

Present neutrino oscillation data fix the mass splittings among the three neutrinos, leaving 
arbitrary the overall scale of neutrino masses. The latter could be as large as an electron volt 
or so without conflicting with cosmology, beta decays and neutrinoless double beta decays. 



given current uncertainties in nuclear matrix element calculations 211. However, in the 

n 

BRpV model, neutrino masses are strongly hierarchical p, specially the lightest neutrino 
mass. As a result the possibility of quasi-degenerate neutrinos is not realized in this model 
and, correspondingly, we will not discuss it any further. In what follows, whenever we fix 
the BRpV parameters, we take the tree level neutrino mass as mj^g = 0.05 eV, the current 
atmospheric best fit value from Ref. jisl . and fix the value of the remaining BRpV 
parameter at representative values 0.22 GeV and 7 x 10~^ GeV. 
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The presence of BRpV induces a mixing between the neutrinos and neutrahnos, giving 
rise to i?-parity violating decays of the LSP. In our model, the lightest neutralino presents 



Xi Jyj^i^ [Ul- The expected Xi lifetime and decay lengths depend both on the magnitude 
of i?-parity breaking parameters and the chosen values of the SUGRA parameters. 

Fig. ^a) shows the lightest neutralino decay length as a function of its mass for Aq = 0, 
yU > 0, tan/5 = 3, and = 0.22 GeV. In this figure, the solid lines stand for the tree-level 
neutrino mass rriiy.^ = 0.05 eV, corresponding to the best fit atmospheric scale as given in 
Ref. 0, and we chose tuq = 200 and 700 GeV. The bands in these figures were obtained by 
taking m^^ within the 3a allowed atmospheric mass splitting of Ref. jla]. Fig. ^b) presents 
the LSP decay length as a function of the neutrino mass, for different lightest neutralino 
masses and the parameters used in Fig. ^a). 

From Figs. ^ we can see clearly that the LSP decay length is shorter for larger neutrino 
and LSP masses, as expected. Furthermore, irrespective of the smallness of the neutrino 
mass indicated by the atmospheric oscillation data, the LSP (x?) decays inside the detector 
in a large portion of the parameter space, specially for small mo. It is interesting to notice 
that the kinks in Fig. ^ are associated to the opening of the LSP decays into on-shell gauge 
bosons. 

Figure |21 contains the Xi branching ratios as a function of mi/2 for e,- = 0.22 GeV, 
rriu^ = 0.05 eV, Aq = 0, and yU > 0. As we can see from this figure, the decay Xi ~^ ^bb 
dominates at small tuq and large tan (3. This decay channel arises from an effective coupling 
A333, induced by the neutrino-neutralino mixing, which is proportional both to the bottom 
Yukawa coupling and to the bilinear parameter er in our model. Therefore, this decay is 
enhanced at small mo due to the lightness of the scalars that mediate it and at large tan f3 
due to enhanced Yukawa couplings. Both features are apparent in the panels of this figure. 
Moreover, whenever Xi ~^ '^bb is not the leading channel, the LSP decays mostly to rud, 
with a sizable branching fraction Xi ~^ ruii. 

In Fig. ini we present the dependence of the Xi branching ratios with respect to the im- 
parity violating parameter e,-. The importance of the ubb decay mode increases for large e,-, 
since the effective coupling A333 is proportional to er- Moreover, Fig. |3] shows that the LSP 
branching ratio into ubb decreases with increasing mj^g in the parameter regions where this 
mode is sizeable. 



leptonic deca^ 




semi-leptonic ones Xi ~^ ^QQ or ^QQ^ ^^'^ the invisible mode 
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One comment is in order here, concerning the naturalness of the above branching ratio 
predictions. This issue depends somewhat on the assumptions about the supersymmetry 
soft breaking terms. We are tacitly assuming them to be universal at some unification scale. 
This is the usual practice and is adopted here to ensure simple comparison with the SUGRA 
i?-parity conserving results for the trilepton signal obtained in Ref . . In the BRpV model 
universality plays another important role, namely, it ensures "calculability" of the neutrino 
mass mj,3, by the renormalization group evolution. One finds in this case that the smallness 
of m^^g required by experiment follows naturally from the approximate "alignment" of BRpV 



parameters discussed, for example, in Ref. jSlllSl- In such scenario one has that the short 
LSP decay path is indeed technically natural, despite the small neutrino mass. However 
one should stress that strictly universal boundary conditions are not at all essential for the 
consistency of the results presented here. 

Summing up, we can say that in the BRpV model the Xi decays mainly into rud for large 
mo or small er, while its decays are dominated by ubb at small mo and large e,- (tan/3). As 
we will see in what follows this has important implications for the trilepton signal. 

III. SIGNAL, BACKGROUNDS, AND SELECTION CUTS 

In i?-parity conserving scenarios, trilepton production at the Tevatron proceeds via pp — >■ 
X2xf with xt ^^Xii X2 ~^ ^^^~Xii ^-nd the LSP (x?) leaving the detector undetected. 
The main SM backgrounds for trilepton production are displayed in Table IH In order to 
suppress these backgrounds, we have imposed the soft cuts SC2 defined in Ref. Il3|, which 
were tailored for scenarios containing soft signal leptons coming from r decays 23[: 

CI: We required the presence of three isolated leptons (e or /i) with a hadronic Et smaller 
that 2 GeV in a cone of size Ai? = 0.4 around the lepton |^|; 

C2: We required the most energetic lepton to satisfy |?7^J < 1.0 and Pt(^i) > H GeV; 

C3: The second (third) most energetic lepton must satisfy l^^^j/sl < 2.0 and (^2(3)) > 7 
(5) GeV; 

C4: We required the missing transverse energy to be larger than 25 GeV; 
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C5: We vetoed events exhibiting a pair with an invariant mass smaller than 20 GeV 
and larger than 81 GeV (this avoids both Z boson and QED contributions); 

C6: We vetoed events with a transverse mass of a charged lepton and missing transverse 
energy between 60 GeV and 85 GeV in order to suppress W decays. 



In our analysis, the signal and backgrounds were generated using PYTHIA |2^, except for 
the WZ*{-y*) which was computed using the complete matrix elements 3]- Furthermore, 
we also simulate experimental resolutions by smearing the energies, but not directions, of 
all final state particles with a Gaussian error given by AE/E = 0.7/y/E [E in GeV) for 
hadrons and AE/E = 0.15/v^ for charged leptons. 

The cross section for the SM backgrounds after cuts are shown in Table HI Notice that 
.e.,t= d.fe fro. o... . Ref. Q became ..e .ad„ni.at.n p.oced.e. 

used in PYTHIA and IS A JET are different. Moreover, we can see that the most important 
background is the ti production, accounting for 70% of the total background. We further 
tested our code by verifying that our results for the i?-parity conserving signal agree with 
the ones in Ref. [13|. 

In our BRpV model there are more SUSY reactions that can contribute to the trilepton 
signal than in the i?-parity conserving case, since the Xi decays can give rise to charged 
leptons. In the parameter range interesting for neutrino physics i?-parity violating decay 
modes are only important for the LSP decays and can safely be neglected for the other SUSY 
particles. We have incorporated all possible decay modes of the lightest neutralino and the 
second lightest neutralino taking fully into account large tan(/3) effects in PYTHIA, leaving 
the other decay modes unchanged. Assuming that gluinos and squarks are too heavy to be 
produced at the Tevatron, we considered the following processes: 

pp-^it , H , x-x;(^(j) = l,2) , xtXi ,and x-X^(^ = l,2). 

The Xi decays can contain charged leptons, and therefore, we should also analyse multi- 
lepton (> 4£) production. In order to extract this signal, we applied the cuts CI, C3, C5, 
and C2 but accepting leptons with \rj{(!^)\ < 3. We also required the presence of an additional 
isolated lepton with px > 5 GeV and demanded the missing transverse energy to be larger 
than 20 GeV. The main SM backgrounds for this process are tt, WZ, and ZZ productions 
whose cross sections after cuts are shown in Table iHl 
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We can see from Fig. that the hghtest neutrahno might not decay inside the detector 
depending on the point of the parameter space. If it decays inside the tracking system, it can 
give rise to spectacular events exhibiting displaced vertices without an incoming track asso- 
ciated to it. In our analyses, we did not look for displaced vertices since the corresponding 
backgrounds depend upon details of the detector. This possibly is a conservative hypothesis 
since this kind of events should present a small background, leading to a larger reach of the 
Tevatron. Nevertheless, we kept track of the position of the neutralino decay and accepted 
events where the neutralino decays inside the tracking system, rejecting all events where 
one of the neutralinos decays outside a cylinder around the beam pipe of radius 0.5 m and 
length 2.0 m; we named this requirement C7. Again, this is another conservative estimate 
of the expected BRpV signal. 

We investigated the regions of the rriQ (g) mi/2 plane where the trilepton and multilepton 
signals can be established at the Tevatron for integrated luminosities of 2 fb~^ and 25 fb~^ 
and fixed values of Aq (=0), tan/5, sign(yu) (> 0), er, and m^g (= 0.05 eV). We evolved the 
renormalization group equations starting at the unification scale (few xlO^^ GeV) for the 
soft parameters, rejecting the points where either the electroweak symmetry is not correctly 
broken, or which exhibit particles with mass excluded by present experimental data jj22|. In 
our analysis, we have employed the Poisson statistics, except when the expected number of 
signal events is large enough to justify the use of the Gaussian distribution j^]. We exhibit 
our results in the mo(^mi/2 plane, denoting by black circles the theoretically excluded points, 
and by white circles the experimentally excluded by sparticle and Higgs boson searches at 
LEP2 I22I. The black squares represent points accessible to Tevatron experiments at 5a 
level with 2 fb^^ of integrated luminosity, while the white squares are accessible with 25 
fb^^. Points denoted by diamonds are accessible only at the 3a level with 25 fb~^, while the 
stars correspond to the region not accessible to the Tevatron. 

The trilepton cross section is always dominated by the XiXi ^^'i X2X1 productions, 
with the second process being about 20% to 50% larger than the first. For example, these 
reactions are responsible for approximately 99% of the cross section at large mo. Note that 
the first process can contribute to the trilepton signal only in i?-parity violating scenarios. 
Moreover, at small and moderate tuq 400 GeV), the xlxf production is responsible for 
approximately 5-10% of the signal cross section and the production of sleptons also gives a 
sizable contribution. 
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To study the multilepton signal we choose three representative parameter regions. The 
best scenario corresponds to the case where tan (3 and e,- are small. We subsequently relax 
this optimistic assumption by considering separately the cases where either (but not both, 
since in this case the signal is too small) or tan/5 are large. 

In Fig. El we present the region of the plane that can be probed at the Tevatron 

for tan/? = 3 and BRpV parameters e,- = 7 x 10~^ GeV and rrii,^ = 0.05 eV. For these values 
of the parameters, the signal cross section is dominated by xtX2 production followed by 
xt ^ Xi^^'^j X2 ~^ XiQQ S'lid the LSP Xi decays mainly into rud. As expected, the long 
lifetime of the neutralino reduces considerably the signal in the shaded area of Fig. after 
we apply the cut C7. Therefore, we might be able to further probe this region by looking 
for displaced vertices, and consequently enhance the Tevatron reach. 

It is interesting to compare our results presented in Fig. Elwith the ones in Ref. First 
of all, the presence of BRpV interactions reduces the Tevatron reach in the trilepton channel 
for small values of rriQ. This happens because there are some competing effects in this region 
of parameters: on the one hand there are new contributions to the trilepton process due to 
LSP decay and on the other hand, the decay of the neutralinos produce a larger hadronic 
activity, worsening the lepton isolation, and reducing the missing Et compared with the 
MSSM case. Besides that, the leptons from the Xi decay can give rise to additional isolated 
leptons which can contribute to the trilepton signal or, alternatively, can suppress it due to 
the presence of more than three isolated leptons. The last effect and the larger hadronic 
activity reduce the trilepton signal at small mo in the BRpV model. In contrast, as can be 
seen from Fig. El the trilepton reach at large mo always tends to increase with respect to the 
MSSM expectation. This follows from the fact that the drastic reduction of the X2 branching 
ratio into leptons at large mo is compensated by the additional production of charged leptons 
in Xi decays. Since these extra leptons come from tau decays, it is important to adopt cuts 
which increase the acceptance of soft leptons. The cuts we used satisfy this requirement. 

In Fig. in^a) we present the Tevatron reach in the multilepton (four leptons or more) 
channel for the same parameters adopted in Fig. El As we can see, the multilepton reach is 
larger than the trilepton one, increasing the discovery potential at large values of mi/2- For 
instance, the Tevatron reach at large mo is ~ 225 GeV in our BRpV model while it is of the 
order of 150 GeV in the MSSM. Moreover, unlike the trilepton signal, the discovery potential 
at small mo is also increased with respect to that of the MSSM. In this region it is clear 
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that the reduction of the trilepton signal is largely due to the presence of additional isolated 
leptons. As in Fig. the shaded area represent the region where displaced vertices could be 
used to further increase the sensitivity to BRpV. In Fig.lHfb) we present the combined results 
for the trilepton and multilepton searches using the chi-square criteria. It is interesting to 
notice that the presence of -R-parity violating interactions leads to a 5cr SUSY discovery even 
at large mo, a region where the usual i?-parity conserving SUGRA model has no discovery 
potential at all. Notice the importance of combining trilepton and multilepton signals to 
achieve this conclusion. 

Let us now consider a second scenario where the Xi decays predominantly into u^bb at 
low rriQ. As seen in Fig. IHl the dominance of this decay channel happens for large e,- values 
since the LSP decay is scalar mediated. The z/366 decay mode spoils lepton isolation, without 
producing any additional charged leptons and may be regarded, in a sense, as a worse case 
scenario in comparison with the case of small e,- and mg. In order to illustrate this case we 
fixed tr = 0.22 GeV and the remaining parameters as before; Aq = 0, tan/? = 3, fi > 0, and 
= 0.05 eV. From Fig. Owe can see that the trilepton reach is indeed further reduced at 
small and medium mo, as expected. The same happens for the multilepton signal; see Fig. 
IS] Nevertheless, combining these signals still leads to a reach that is better than the MSSM 
one for all mo and mi/2 values. 

Finally, we consider the case of small and large tan/3, say = 7 x 10~^ GeV and 
tan (3 = 35. Fig. IHl displays the Tevatron reach in the trilepton channel for this case, keeping 
the remaining parameters as before; Aq = 0, fi > 0, and m^^ = 0.05 eV. For these parameters, 
the main Xi decay mode is rud. However, there is a sizable contribution of the ubb channel 
at small mo. As expected, the SUSY reach decreases at small mo, specially as we increase 
tan/3. In contrast, there is a slight gain for mo ^ 200 GeV; see Fig. El The similarity 
between the results in Fig. El and Fig. El at large mo can be ascribed to the rud dominance 
of the LSP decay; see Fig. El^b) and Fig. El^d). In contrast the situation is more complicated 
at smaller mo as seen in Fig. and Fig. IHl 

For the last choice of parameters, the Tevatron discovery potential in the multilepton 
channel is larger in our BRpV model than in the MSSM except at small mo, however, it 
is similar to the first case we analyzed; see Fig. iTUT a) and Fig. IHl^a). This feature survives 
the combination of three and multilepton signals, as can be seen from Fig. Efb) and Fig. 
irUT b). In contrast, the combined analysis does make a difference in the small mo region. In 



11 



all cases, however, the reach of the BRpV model is larger than the MSSM one. 

IV. CONCLUSION 

We have analyzed the production of multileptons ( > M with £ = e or yu) in the simplest 
supergravity model with violation of R parity at the Fermilab Tevatron. In this model, an 
effective bilinear term in the superpotential parameterizes the explicit breaking of R parity. 
Despite the small i?-parity violating couplings needed to generate the neutrino masses indi- 
cated by current atmospheric neutrino data, the hghtest supersymmetric particle is unstable 
and can decay inside the detector. This leads to a phenomenology quite distinct from that 
of the /2-parity conserving scenario. We have quantified by how much the supersymmetric 
multilepton signals differ from the i?-parity conserving expectations, displaying our results 
in the mo ® Tni/2 plane. We have shown that the presence of bilinear i?-parity violating 
interactions enhances the supersymmetric multilepton signals over most of the parameter 
space, specially at moderate and large mo. These topologies are useful not only for discovery, 
but also to verify whether R parity is conserved or not. 

Adopting the hadronization procedures used in PYTHIA, we have first reproduced the 
results for the trilepton signal expected in the conventional i?-parity conserving supergravity 
model. We have found good agreement with the results of Ref. which adopts the ISA JET 
event generator. We have then shown how the presence of BRpV interactions leads to a 
small suppression of the trilepton signal at small values of mo irrespective of the value of 
BRpV parameter e^. This is due to the x? decay into vhb. However, the x? decays lead to a 
drastically extended reach at large mo as a result of the LSP decay into rud. Moreover, the 
presence of additional isolated leptons in the signal allows us to look for multilepton events. 

We have demonstrated that combining the trilepton and multilepton searches increases 
the Tevatron Run II sensitivity for most of SUGRA and -R-parity breaking parameters. 
Note, however, that for neutrino masses in the range indicated by current atmospheric data 
one has a gap between the x? niasses that can be probed at LEP2 (up to 40 GeV or so) 
and those that can be studied at the Tevatron (above 70 GeV or so): within this range the 
Xi decay length is rather large, requiring the study of other topologies, like the presence of 
displaced vertices in the tracking system. It is interesting to notice that we can search for 
SUSY signals also in the low mo region by looking for events exhibiting multi jets + lepton 



+ missing transverse momentum 
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In the present paper we have confined ourselves to the case in which the hghtest neutrahno 
is also the lightest supersymmetric particle j^^], the most likely possibility if we adopt the 
simplest set of supersymmetry soft breaking terms, universal at some unification scale. We 
also have focused on the case where we have only one generation and this is chosen to be the 
third. This is done first for simplicity. Second we adopt this choice as a worse-case scenario. 
In other words, in those parameter regions where our multi-lepton signal can be discovered in 
the present one-generation approximation, the inclusion of additional generations can only 
improve our result. In contrast, in those regions where our results are negative, the situation 
is totally inconclusive in the sense that a full fiedged analysis including all generations might 
reveal that the signal can also be detected in part of those regions. Therefore our results are 
robust, in the sense that the inclusion of additional generations would imply new sources 
of leptons, specially muons. The analysis is substantially more involved, however, than the 
one considered here and will be taken up elsewhere. 
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BG (fb) 


a(fb) 


WZ (Z ^ tt) 


0.17 


W*Z*,W*7 ^ III 


0.12 


W*Z*,W*7 ^ ZZ'Z^ 


0.15 




1.15 


2*2* 


0.05 


total 


1.64 



TABLE I: Background cross sections in fb for the trilepton signal at the Tevatron Run II after 
kinematical cuts discussed in the text. 



BG (fb) 


a (fb) 


WZ 


0.01 


2*2* 


0.10 


tt 


0.16 


total 


0.27 



TABLE II: Background cross sections in fb for the multilepton signal at the Fermilab Tevatron 
Run II after kinematical cuts discussed in the text. 
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FIG. 1: Xi decay length versus LSP mass for Aq = 0, fi > and tan/? = 3 for (a) fixed BRpV 
parameters: = 0.22 GeV, m^^ = 0.05 eV (solid lines) and current atmospheric 3a band (shaded 
bands); (b) as a function of the neutrino mass m^^ for tuq = 100 GeV and several values of Xi 
masses . 
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FIG. 2: Xi branching ratios as a function of mi/2 for = 0, ^ > 0, e,- = 0.22 GeV, and = 0.05 
eV. The sohd hnes denote Xi i^s^b (squares); Xi — ^ Tud (circles); and Xi — ^ fs^"*"^" (stars). The 
dashed hnes denote Xi — ^ invisible (squares); Xi (circles); and Xi — > ^^399 (stars) . 
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(a) tan|3 = 3 Mj = 50GeV (b) tan|3 = 3 Mo = 800GeV 




e e 

T X 

FIG. 3: Xi branching ratios as a function of e,- for j4o = 0, ^ > 0, and 771,^3 = 0.05 eV. We fixed 
"^1/2 = GeV in the case of tan/3 = 3 and mi/2 = 125 GeV for tan/3 = 35. The hnes are as in 

Fig.m 
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FIG. 4: Xi branching ratios as a function of niy^ for = 0, ^ > 0, and e,- = 7 x 10~^ GeV. We 
fixed ?7ii/2 = 175 GeV in the case of tan/? = 3 and = 125 GeV when tan/3 = 35. The Unes 
are as in Fig. [2 
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FIG. 5: Reach of Fermilab Tevatron Run II using the trilepton signal in the mo iX" w,i/2 plane 
for A[) = 0, tan (3 = 3, fi > 0, er = 7 x 10~^ GeV, and rriu^ = 0.05 eV. The black circles are 
theoretically excluded, while the white circles are experimentally excluded by sparticle and Higgs 
boson searches at LEP2. The black squares denote points accessible to Tevatron experiments at 
5(7 level with 2 fb~^ of integrated luminosity, while the white squares are accessible with 25 fb~^. 
Points denoted by diamonds are accessible at the 3a level with 25 fb~^, while the stars correspond 
to the region not accessible to Tevatron. The long lifetime of the neutralino reduces considerably 
the signal in the shaded area, however, it suggests that the sensitivity can be improved by looking 
for displaced vertices. 
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FIG. 6: (a) Reach of Fermilab Tevatron Run II in the 4 or more lepton channel, (b) Combined 
trilepton and multilepton results. All parameters and conventions were chosen as in Fig. IHI 
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FIG. 7: Reach of Fermilab Tevatron Run II using the trilepton signal in the mo rni/2 plane for 
^0 = 0) tan/3 = 3, ^ > 0, e,- = 0.22 GeV, and 171^3 = 0.05 eV. The conventions are the ones in Fig. 
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FIG. 8: (a) Reach of Fermilab Tevatron Run II in the 4 or more lepton channel, (b) Combined 
trilepton and multilepton results. All parameters were chosen as in Fig. d 
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FIG. 9: Reach of Fermilab Tevatron Run II using the trilepton signal in the mo rni/2 plane for 
^0 = 0) tan f3 = 35, ^ > 0, er = 7 X 10~^ GeV, and m^^ = 0.05 eV. The conventions are as in Fig. 
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FIG. 10: (a) Reach of Fermilab Tevatron Run II in the 4 or more lepton channel, (b) Combined 
trilepton and multilepton results. All parameters were chosen as in Fig. IHl 



